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1. Introduction 

Coherent transports of electrons through a quantum dot (QD) embedded in an Aharonov- 
Bohm (AB) ring have been widely studied in recent years. For such systems, interference be- 
tween the localized states in the QD and the continuous states in the direct scattering channel 
causes asymmetric resonance peaks in the conductance/'^ the so-called Fano resonances.'^ On 
the other hand, strong Coulomb interaction induces the Kondo effect, i.e., the formation of a 
singlet bound state among localized and conduction electrons. In case of transport through 
single QD, it is seen as a plateau in the resonance tunneling peak, where the conductance 
reaches the unitary limit value (2e^//i). Thus, when the Coulomb interaction is strong enough, 
an interplay between the Fano and Kondo effects is expected in the AB ring with a QD. 

In the linear response regime, the interplay was theoretically studied by Hofstetter et 
ai, employing the numerical renormalization group (NRG) technique,'^ and they predicted 
suppression of the Kondo plateau by the Fano resonance (the Fano-Kondo resonance). Such 
a behavior has been observed by recent experiments.^'^ They also studied the AB oscillations 
in the plateau regime"^ and showed that the phase of the AB-oscillation is almost fixed at 
7r/2^ and its frequency is almost doubled, both relative to the AB-oscillation at higher gate 
voltage. 

Outside the linear response regime, the conductance maximum at zero bias voltage^' ^ 
characterizes the Kondo effect in the transport through a single QD. The conductance drop 
at finite bias voltage is considered to be due to inelastic scattering involving electrons with 
energy between the chemical potentials of the two electrodes in case of infinite Coulomb 
repulsion. An AB ring with a QD is predicted to have a conductance minimum or maximum 
at zero bias voltage depending on the magnetic field. In spite of these researches, we think 
that the nonequilibrium properties of an AB ring with a QD particularly in case of large but 
finite Coulomb interaction are not fully understood partly because of the lack of complete 
characterization of nonequilibrium ensembles. 

On the other hand, the method of C*-algebra has been widely applied to rigorous investi- 
gations of nonequilibrium steady states (NESSs) of infinitely extended quantum systems. 
The method was originally developed to deal with equilibrium statistical mechanics^^"^^ 
and, then, found to be useful also for studying nonequilibrium properties. ^^"^^ Particularly, 
Frohlieh, Merkli and Ueltschi^"*^ have shown the existence of NESS and its characterization for 
junction systems, which consist of several free-fermionie reservoirs mutually interacting via 
bounded interactions. Their results cover the AB ring with a QD in case that the electrodes 
are free fermion systems, but explicit characterizations of NESS are not given. 

In this paper, we investigate the transports in an AB ring with a single-level QD, em- 
ploying the finitc-f/ slave-boson mean-field (SBMF) approach of Kotliar and Ruckcnstein^^ 
extended to nonequilibrium steady states. The NESS with respect to the mean-field Hamilto- 
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nian is constructed based on the C*-algcbraic approach as in the noninteracting case.^^' The 
infinite-C/^^"^^ and finite- ^7^''"^^ SBMF approaches in combination with the Keldysh Green 
function method are successfully applied to the study of strongly correlated systems. In par- 
ticular, the main features of the Kondo effect in the transport through a single QD such as 
the conductance plateau and the zero-bias maximum have been well described by the finite-C/ 
SBMF approach.^*^'"*^ In contrast to these works, we focus on the finite-C/ case and treat the 
mean-field NESS in terms of the C*-algebraic approach, which leads to a global characteri- 
zation of NESSs and, we believe, has an advantage compared to the Keldysh method in this 
respect. 

The rest of this paper is arranged as follows. In Sec. 2, after describing the model of the 
AB ring with a QD and reviewing the slave-boson description of Kotliar and Ruckenstein, the 
nonequilibrium SBMF approximation is explained. Relations among mean-fields and Lagrange 
multipliers are derived from the equations of motion of the auxiliary bosons within Dirac's 
formulation of constrained dynamics in contrast to the original SBMF approach, which 
is based on the saddle-point approximation of the path integral representation of the free 
energy. The self-consistent equations are closed by evaluating the averages of certain fermionic 
variables with respect to a NESS for the mean-field Hamiltonian, which is constructed with 
the aid of the C*-algebraic approach. Section 3 is devoted to the summary of the self-consistent 
equations and the derivation of the average current. Within this approximation, the Coulomb 
interaction is shown to result in the renormalization of the dot energy and the reduction of the 
dot level-width. In Sec. 4, we study the condiictancc in the linear response regime and show 
that (i) the Fano-Kondo resonances and AB oscillations obtained from the SBMF approach 
are in good agreement with the results from the NRG technique by employing twice large 
Coulomb energy.^ In addition, the SBMF approach is shown to give a simple view to the 
underlying physical processes. In Sec. 5, we study the differential conductance at fixed average 
number of dot electrons. We have found that (ii) the Fano-Kondo resonances arc deformed 
by the bias voltage so that the resonance peaks have the tendency to split into two peaks and 
that (iii) the differential conductance may have zero-bias maximum or minimum depending 
on the transmission via the direct tunneling between the electrodes. The former is due to the 
narrowing of the dot level-width by the Kondo effect and, at low bias voltage, can be observed 
as an increase of the conductance plateau, which is consistent with the recent observation by 
Katsumoto et. al. for the T-shaped QD systems.^ Section 6 is devoted to concluding remarks. 
In Appendix, we explain the details of the construction of asymptotic field operators necessary 
for the construction of NESS. 
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2. Nonequilibrium Slave-Boson Mean-Field Approach 

2.1 Model and slave-boson approach 

Our system is described by the Anderson single impurity Hamiltonian^^ with the direct 
channel between left and right electrodes: 

H = Hl + Hr + Hd + Ht + Hlr. (1) 

The first and second term are Hamiltonians of the two-dimensional electrodes 

Hr = ^ dkLOkralar^kar (r = L or R) (2) 
a 

where ujkj. is the single-electron energy: = ^kR — eV = h'^k'^/{2m) — eV/2 with V the 
bias voltage and m the effective mass, and al^^^akar) is the creation (annihilation) operator 
of the electrode electron with wave number k and spin a. The third term Hd describes the 
single-level dot with the Coulomb interaction U : 

Hd = ^ eaclca + Un'^ui (3) 

a 

where Ca{c^) is the creation (annihilation) operator of electrons in the dot, = c^c^ the 
number operator and the single-level energy. The tunneling between the dot and electrodes 
is given by 

Ht = ^ dk[ukLa\^i^c^ + UkRa\^j^c„ + (H.c.)] (4) 

and the transmission through the direct channel by 

Hlr = Y.jj dkdq[WukLUgRe''^al^j^agaR + (H.c.)] , (5) 

where the tunneling matrix elements UkL, UkR and WukLUgR (hence W) are taken to be real. 
Because the resonance shapes do not strongly depend on the functional form of the tunneling 
matrix elements, one can assume that the interaction Hamiltonian H^r is separable as above 
without loss of generality. The AB phase (p is related to the magnetic flux # penetrating the 
ring via ip = e^/{hc) with c the velocity of light. 

In order to examine the transport through the strongly correlated QD, wc employ the 
finite-?7 slave-boson mean-field (SBMF) approach proposed by Kotliar and Ruckcnstcin.^^ 
In order to describe the doubly occupied state as a bosonic single particle excitation, they 
have introduced four auxiliary bosonic fields e,Po- {(^ =T)i) and d, respectively, labeling the 
empty, singly occupied and doubly occupied states, and a fermionic field f„, describing the 
fermionic feature of the QD electrons. The expression of the Hamiltonian (1) in terms of the 
new variables is obtained by replacing the original QD electron operators as follows: 

ctrCfj > fa fa 

n^ni — Sd 
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where is a bosonic variable defined by 

= (1 - S}d-v\v,)-^l\e^v,^v\m - e^e -4p^)-V2. (6) 

The additional degrees of freedom simplify the Coulomb interaction IJn'^ni as C/d^d, but, at 
the same time, introduce vectors not describing physical states. Then, physical vectors l^p/i) 
are picked up by a set of constraints 

HpIp. + d)d-flu)\%h) = 

H%h) ={e^e + dW + J^PlPa - l)l^^>/^) = (7) 

cr 

As easily seen, the subspace of the enlarged Hilbert space defined by (7) is equivalent to the 
original Hilbert space. Note that, since (f)ia,4>2 and H commute with each other, no further 
constraints are necessary. 

As discussed in detail by Dirac,^^"^^ corresponding to the constraints (f)ia\^ph) = 
and 4>2\^ph) = 0) Q'-number Lagrange multipliers A^^ and A^^^ should be introduced 
to the Heisenberg equation of motion and, with respect to the effective Hamiltonian 

any dynamical variable A satisfies the standard equation of 

motion as a state equation 

rIA -i ~ 

— \^^^) = -[H,A]\^,^) (8) 
where \^ph) is an arbitrary physical state and the explicit expression of H is 

H =Hl +Hr + Hlr + Y. I dk[zl{UkLflakaL + UkRfWaR) + {H.C.)] + USd 

+ J2 '-f-f- + E ^ (p-p- + - + ^^'^ (^^^ + + E p-p- - 1) • (9) 

Note that, under the constraints (7), the equation of motion (8) is equivalent to the Heisenberg 
equation of motion with respect to the original Hamiltonian (1). 
2.2 Nonequilibrium mean field approximation 

The key idea of the SBMF approximation is to replace all the slave-boson operators and 
the Lagrange multipliers by their average values, which will be denoted as e,pa,d, and 
A(2), and to impose, instead of (7), weaker constraints among averages: 

|ep + |d|2 + ^|p,|2 = i (10) 

a 

{f^Ja)ss = \df + \Pa\\ (11) 

where {faf(j)ss stands for an average of fafa- This approximation can be extended to nonequi- 
librium steady states (NESSs) by regarding e,Po-5 c?) \ A^^^ and {fafa)ss as the NESS aver- 
ages with respect to the mean field effective Hamiltonian: 

H^ff = H, + U\d\' + ^ A«(|p.p + \d\^) + A(2) (|e|^ + \d\' + ^ \p^\' - l) (12) 
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where Hp is the fermionic part: 
Hf = Hl + Hr + HlR + J2 f dk[z;{UkLfWaL + UkRfW.R) + {H.C.)] + ~^-flU , (13) 

= — is the effective energy of the QD level and the mean value z„ is obtained from 
(6) by replacing the bosonic operators to their averages. The characterization and properties 
of NESS with respect to i^e// will be discussed in the next subsection. 

We have seven averages to be determined and the constraints provide three conditions. 
Hence, four more conditions are necessary. Originally, Kotliar and Ruckenstein^^ derived them 
by the saddle-point approximation for the equilibrium free energy, but this approach can not 
be applied to NESS. Here, we derive those conditions from the equations of motion of the 
slave boson fields. For example, d obeys the equation of motion: 
d f 

^^dt^ = - X] / dk{[zl,d]{ukLfWTL + UkRflakaR) + [Z(T, A{'^kLa\^LU + UkRoiaRU)} 

+ (C/ + XAW+A(2))d. (14) 

a 

By replacing the slave-boson operators and Lagrange multipliers to their mean values and by 
evaluating NESS averages of the fermionic operators, one gets 

where Mo- = J dk{ukLfaakaL + UkRfaakaR)ss- Here, the average of [zcr,d] is evaluated as 
[za,d] = [(1 - d^d-plp^r'^^d] {e^pa+pU)il - eh-plp-,)-^'^ 

= E %n)T' i(d^d + plp,r,d] (etpo+4d)(l-ete-4p<,)-V2 

n=l ^ 

= - E (2J.\' T.(d^d + plp,r-'-'d{SUplpa)\e^Pa+pW^ - e^e-pU,)-'/' 



n=l 



dz: 



(2n)!! ^' ' ' ' ' ' ' ' ' dd* 



and that of [zl, d] as —-^r-^- Similarly, from the equations of motion of e and p^-, one has 

dd* 

E(f§^^ + |^K)+^^'^^ = o (16) 

E (||^^ + '^^^r) + (A(^) + A(^))P. = (17) 

When one neglects the Zeeman effect as in this article, a straightforward investigation 
shows that e,Pa,d and Za can be taken as real positive and, then, Mg- is real. Therefore, by 
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setting d = VD, Pa = VPa, g = ^/E, the constraints (10), (11) and the conditions (15), (16), 
(17) cast into self-consistent equations for D and A^': 



^{l-D-P^){l-E-P^) 
2^5^M, + C/ = (20) 



d d d d 
where ^ = ^ + ^ ~ XI ^P" ^"^^ ^"^^ should set = {flU)ss -D,E = D + l- 

Ylaif^Mss after evaluating all the derivatives. 

2. 3 Nonequilibrium steady states for the mean field Hamiltonian 

Now we turn to the averages of fermionic operators which are evaluated at a nonequihb- 
rium steady state (NESS) with respect to the fermionic part Hp of the effective Hamiltonian. 
Since we are interested in transports through mesoscopic contacts, we consider a NESS such 
that, very far from the contact, each electrode recovers equihbrium with individual tem- 
peratures and chemical potentials. According to the C*-aIgebraic approaches,^^'^^'^^^^'^^ an 
average of an observable A with respect to such a NESS is can be evaluated as 

{A)ss = Mm {e'^'Ae-'^\ (21) 

where {A)q = Tr{piA) is the average with respect to a local equilibrium state pf.^'^ 

= }Lg-PL{HL-IJ.LNL)g-PR{Hi{-ij.i{Nii) _ ^22) 

In the above, S is the normalization constant, /S^, Mr and A^^ (^ = ^^P) are, respec- 
tively, the inverse temperature, chemical potential and the total number of particles: 
Nr = J2a I '^kar^kar of each electrode. 

The existence and characterization of the NESS averages can be heuristically explained in 
terms of the M0ner operator: = lim g'^HFt/h^-iHot/h -y^rhere Hq = Hi + Hr. Indeed, since 

t— oo 

the initial density matrix pi commutes with Hq, one has 

{A)ss = Hm l^e-'^ot/h^iHpt/h^^-iHFt/h^iHot/h~^^ ^ ^j^^ (^^iHot/h^-iHFt/h^^iHpt/h^-iHot/h^^ _ 
t— »+oo t— »— oo 

This seems to suggest {A)ss = {^^ACI)q, but, as in the derivation of the Gell-Mann Low 
relation,^^ the change of normahzation should be taken into account and one obtains 

{AUs = . (23) 

According to the rigorous approaches,^^'^^'^^-*^^'^^ the existence of {A)ss for an aribtrary finite 
observable A requires the absence of bound states of Hp, which is assumed throughout this 
article. In this case, = 1 and 11 = Cl^^ is a projection onto the ortho-complement of the 
states involving dot electrons, namely, one should have 11 = (1 — /|/|)(1 — fjfi)-^'^ 
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Simple characterization of the NESS is given in terms of the asymptotic fields. 
According to the scattering theory,^^ the M0ller operator Q, generates an 'in' field: 
pl. = lim e'^Ft/n^-iHot/hJ 'iHot/h^-iHFt/h ^ nai.n^, which creates 'in' states from 

t — '' — CO 

the vacuum and satisfies the same anti-commutation relations as afeo-r-. Indeed, one has e.g., 
= S{k — k')dcra'Srr'^^^^^^ = S{k — k')5cr(j'Srr' ■ 

Also the intertwining relation Hp^ = CIHq shows that they are normal modes of Hp'- 
[HF,Pkar] = ^[Ho,akcrr]^^ = —^^kr^cikcrr^^ = —^krPkar- The NESS averages of a product of 
'in' fields are then calculated as follows: 

il^kiairi(^k2o-2r2/~^^s'^'-i'''s^^t'^i'''i^ -^^ (^ ^ki(Tirif^k2(T2r2^^'3'^3r:i ^)o/ (n)o 

= (n4iairin42<72r2nafc3^3r3nafc4<,4^4n)o/(n)o 

= (4i(7iri'^l2(T2r2"'fci'^3r3afc4(T4r4n^)o/ (n)o 

~ (^Ilo-iri'^l2CT2r2^'^'3'^3'-3'^fc40-4r4)o(n)o/(n)o = ('^IiCTiri'^l2CT2r2^'^30'3r30'A;40-4r4)o ; 

where we have used [afccrr^n] = 0, = 11 and the fact that akar and fa (hence H) are 
independent in the initial state pi. Repeating the same arguments and reminding that the 
initial state pi satisfies Wick's theorem for Ukar- one finds that the NESS satisfies Wick's 
theorem for the 'in' fields P^ar and, in particular, the two-point functions are given by 

{Plar0k'a'r')ss = Fr{uJkr)S{k - k')S^^>Srr'- (24) 

where Fr{e) = l/(^e^r{e-ftr) _|_ stands for the Fermi distribution function of each electrode. 
Reminding that the 'in' states are free particle states at t = — oo, eq.(24) simply implies 
that the particles originally coming from the left (right) electrode carries its initial statistical 
information such as the temperature and chemical potential. 

In order to close the self-consistent equations, one task is left, namely the construc- 
tion of the 'in' fields. As shown in Appendix, they are obtained from an operator equation 
[HF,Pkar] = —^krPkar with the boundary condition: akarit) ew{i'^krt/fi) — ^ Pkar as t — > — oo. 
The so-obtained relations can be inverted and the original operators are expressed in terms 
of the 'in' fields: 

fT= yZ f dkzaUkr ^'^l^'^,' ^ ^!'^' (25) 



^ Z^UkrUk'rPk'ar Cai'^k'r)* , zlukrUk'fPk' af Hai^k''. 



'^^ J V UJk'r - ^kr + «0 ko- {uJk'r ) * ^k'f - <^kr + ^0 {oJk'r)* ) 

where r = L or R and we abbreviate L = R, R = L. The auxiliary functions are defined as 

ka{x) = iya{x){x - eg) - y^^Qgrjx) - 2W COS ip^aL{,x)Q,aR{x) / Zg 
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= l + n„,.{x)[W^{x-e„)/zl + 2Wcoii^/za] 
K„{x) = \^We"^-{x-l„)lz„ 
f u^,, z^ 

where Q^cri^) = / dk" — and (pi {'-Pr) stands for ip {—(p). 

J X — LOf^i'j. — lO 

By combining (25), (26) with (24), one can calculate {flfa)ss and in (10), (11), (15), 
(16) and (17) (or in (18), (19) and (20)). This completes the derivation of the self-consistent 
equations. We note that, for equilibrium states, the self-consistent equations derived here 
reduce to those obtained by Kotliar and Ruckenstein^^ via a saddle-point approximation. 

3. Basic Equations 

Prom now on, we assume that the wave-number dependence of the tunneling matrix ele- 
ments is weak so that Ukr = Ur {r = L, R) is constant and that the electrode-dot coupling is 
symmetric: 2'ku\pl = 27rv!j^pR = T, where Pl/r is the density of states of the left/right elec- 
trode. Moreover, we neglect the Zeeman energy and limit our attention to the spin degenerate 
case. Then, the mean fields, the average electron number per spin at the dot as well as the 
original dot-level are independent of the spin quantum number: z^ = zi = z, = X^^^ = A, 
if] 11)33 = = n, and e-f = ej^ = eo. And, the self-consistent equations (18), (19) and 

(20) read as 

V(n -D)(l + D-2n) + ^D{n - D) ^^7) 



V^n(l - n) 

M 



\Jn{l - n) 
dD' 



2 n-D l + D-2n D ^ z{l - 2n) 



l + D-2n V n-D \ n-D ^n{l - n) 



(28) 



4^M + ?7 = (29) 



where 

= -7: — 

2tt 



(l + /x) J dejj^{FL{e) + FR{e)) + ^V^smip J dej^{FL{e)- FR{e)) (30) 

M = ^{1 + z'x) I de^i^i^(F^(6) + Fnie)) 

+ ^^sm^ J deil^f^±^(F^(e) - Ffl(e)) (31) 

with |A(e)p = [(l+.T)(e/r-(eo-A)/r)-|- cos p]'^ + z'^ and x = n'^W'^uj^pL u\pR measures 
the strength of the transmission via the direct tunneling. The chemical potentials pLr of the 
electrodes, appeared in the Fermi distribution functions F,.{e) = l/(e^''('^~'^''^ -|- 1) (r = L,R), 
are taken to be hl = M — ey/2, fiR = p, + eV/2, where p is the average chemical potential. 
The probability D of the double occupation at the dot level and the Lagrange multiplier A 
are obtained as solutions of these self-consistent equations. 
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Now wc turn to the average current at the nonequihbrium steady state. The current 
operator from left to right electrodes is expressed by the time derivative of the number operator 
of the left electrode Nl = Yla ^ '^^^\aL^kaL'- 

= ^ XI / dkukLZaa\^Lf'^ + jj dkdqWukLUqRe"^al^j^aqaR - (H.C.) . (32) 
Its NESS average is given by 

{J)ss = J d-kukLZaial^L fa) SS + JJ dkdqWukLUqRe^'^{al^j;^ag^R)ss- (c.c.) 

= -- / deT{e){FL{e)-FR{e)}. (33) 

where the first term is evaluated as {z^a)j^^j^fa)ss = Za{(il.^i^fa)ss and (24) is used. In this 
equation, T{e) is the transmission probability: 



T{e) 



2-y/x / e eo — A 



z^ \r r 



n 2 



+ cos (f 



+ sm^(p} . (34) 



It is interesting to rewrite the transmission probability in the Fano form:^'^^'^^ 

= T,|±^ (35) 

where E = (1 + x){{e - cq + X)/{z'^T) + [y/x/{l + x)]cos (pj and = 4x/(l + x)^ is 
the background transmission probability.^ The Fano parameter q = [(1 — x)cosip + i[l + 
x) sin(p]/2^/x takes into account the interference effect and determines the peak shape of the 
conductance. The conductance shows the Breit-Wigner type resonance if x = (the real part 
of q is oo) and the anti-resonance if x = 1 (the real part of g is 0). Eq. (35) also indicates that 
the SBMF approach takes into account the effect of Coulomb interaction in two ways, namely 
as renormalization of the dot energy by A: eo — > eo — A, and as reduction of the electrode- 
dot interaction by a factor of 2:^: F z^T. The energy renormalization is partly taken into 
account in a simpler Hartree approximation. The reduction of the electrode-dot interaction 
implies the increase of the lifetime of the dot level and is interpreted as the stabilization of 
the dot level due to strong correlation between the conduction and localized electrons, which 
might be identified with precursor of the formation of a singlet quasi-bound state between the 
two and, thus, with development of the Kondo correlation. 

In the following sections, we investigate transport properties at zero temperature, 
= I/Pr = 0, in detail. 

4. Transports in Linear Response Regime 

In this section, we focus on the linear response regime where the conductance is given by 

2e2 
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and the mean fields D and A are determined by (27), (28), (29), (30) and (31) with 
V = 0. Then, we compare the results with those by the numerical renormalization group 
(NRG) technique.'^ 
4-1 Fano-Kondo resonances 

Figure 1 shows the conductance as a function of the original dot level eo for various 
values of the background transmission probability where the phase difference due to the 
magnetic flux, the electrode-dot coupling and the Coulomb energy are taken, respectively, as 
(/? = 0, r/^ = 0.02 and U/ ^ = 0.6 with ^ the average chemical potential. When Tf, = 0.0, the 
conductance reaches the unitary limit G = 2e^ /h and the well-known Kondo plateau is formed. 
With increasing the background transmission probability (cf. the curves for T5 = 0.1, T5 = 0.3), 
the Kondo plateau is modified by the Fano resonance due to the direct channel transport. 
When Tf, = 0.6, the conductance curve has an asymmetric form, typically observed in the 
Fano resonance, with a large plateau between the peak and dip. The plateau is generated by 
the strong Coulomb interaction. Indeed, as shown in Fig. 2, the plateau region is gradually 
formed by increasing the Coulomb interaction. When the transport through the direct channel 
is dominant (i.e., T5 = 1.0), the destructive interference causes the anti-resonance. Such a 
behavior is experimentally observed in the T-shaped QD system.^' "''^ 

We would like to remark that all the curves in Fig.l are in good agreement with those in 
Fig. 2 of Ref.4, that is obtained by the NRG technique. The parameters of the two results are 
almost the same except the strength of the Coulomb interaction [//F ~ 30, which is almost 
twice larger than that of Ref.4. This indicates that the essential features of the conductance 
such as the shapes or peak positions arc well captured by the SBMF approach in spite of 
neglecting the fluctuations of the bosonic fields and Lagrange multipliers. 
4.-2 Conductance plateau 

Figure 2 indicates that strong Coulomb interaction is responsible for the plateau formation 
in the conductance curve. The underlying physical mechanism becomes more transparent by 
studying the average number of electrons per spin occupying the dot: n, and the probability 
of simultaneously finding two electrons at the dot: D. The left side of the upper panel in Fig. 3 
shows the conductance at = 0.6 and the right side shows the average occupation number 
n, and the double occupation probability D. When the dot level is far above the Fermi energy 
{eo/n > 1.0), no electrons are found in the QD. As the dot level lowers, the occupation number 
n increases and reaches 1/2 at eo/fJ, — 0.9. Within the interval 0.5 < co/a* ^ 0.9, one sees 
n ~ 1/2 and D ~ 0, which implies that the dot is singly occupied and the second electron can 
not enter into the dot owing to the strong Coulomb interaction. In this interval, the electronic 
state does not change and, as a result, the conductance shows the wide plateau. When the 
dot level lowers further, the second electron starts to flow into the dot as indicated by the 
increase of the double occupation probability D, and the conductance again starts to form a 
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peak. 

It is interesting to see the associated behavior of the SBMF parameters. The lower panels 
in Fig. 3 show the reduction factor of the dot level- width and the renormalized dot energy 
eo — A. Outside the plateau region 0.5 ^ eo//i ^ 0.9, ~ 1 and eo — A oc eo + (const.). Hence, 
the dot level-width and the renormalized dot energy are the same as their bare values and 
the behavior of the system is essentially equivalent to that of the nonintcracting case. On the 
other hand, within the plateau region 0.5 < eo//fi ^ 0.9, eg — A is almost constant and 
is very small. The former indicates that the Coulomb interaction modifies the efi'ective dot 
energy seen by the reservoir electrons so that the dot level is kept to be singly occupied. The 
latter indicates that the dot level-width is drastically suppressed by the Coulomb repulsion. 
The formation of the plateau is mainly due to the renormalization of the dot energy. As will 
be discussed in the next section, the effects of the reduction of the dot level-width can be seen 
more clearly in nonequilibrium regime. 
Jj..3 Phase dependence of conductance 

Figure 4 illustrates the phase dependence of the Fano-Kondo resonance. As in the nonin- 
tcracting case, with changing the AB phase from (p = to (p = it, the destructive interference 
between the localized and continuous states turns into the constructive interference and vice 
versa. On the other hand, the plateau of the two curves are almost the same. At ip = tt/2, the 
conductance has a symmetric resonance peak and reaches the unitary limit value as in the 
case of Tb = 0.0 and = (cf. Fig.l). 

We now turn to the relation of the AB oscillations with the Fano-Kondo plateau. It was 
shown that the phase of the AB oscillation at one side of the conductance peak is different by 
vr from that at the other side.^^'^^ Also the strong Coulomb interaction is known to fix the 
maximum position of the conductance at (y9 ~ 7r/2 in the middle of the Fano-Kondo plateau.'^ 
Figure 5 shows the conductance as a function of the AB phase p for various occupation 
numbers n of the dot level. Outside the plateau regime (n = 0.2, 0.8), the conductance shows 
a 27r-periodic AB oscillation and dose not reach the unitary-limit value 2e'^ /h. In contrast, 
in the plateau regime {n = 0.45,0.5,0.55), the conductance reaches the unitary limit 2e^ /h 
near p = tt /2 and is similar to a vr-periodic sinusoidal curve. In other words, in the plateau 
regime, the period of the AB oscillation is almost fixed to be vr and the conductance takes its 
maximum near p = tt/2. This behavior is different from the nonintcracting case where the AB 
oscillation gradually changes its character as the increase of n?^ Note that these features of 
the AB oscillation in the interacting case are consistent with the experiments and theoretical 
predictions of the works mentioned above. 

5. Transports in Nonequilibrium Regime 

It is often argued that SBMF approaches are not appropriate for studying the dynamical 
properties of strongly correlated systems since they can not deal with the effect of charge 
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fluctuations. ' We have shown that the finite-^/ SBMF approach describes weh, at least 
qualitatively, the transports of the AB ring with a QD in the linear response regime as 
compared with the NRG approach.^ So we believe that the SBMF approach could catch 
the essential features of the steady-state transports at finite bias voltage V. In this section, 
we study the properties of the differential conductance at fixed occupation number n: 



where the subscript n means that the F-derivative is taken by keeping n constant. Instead of 
the differential conductance at fixed dot level eg, here we consider Gn because of a technical 
difficulty in solving the self-consistent equations. But the difference is small enough for de- 
scribing the essential features except the region of |/iG„/(2e^)| <^ 1. For example, at (/? = 0, 
T/fi = 0.02, U/fi = 0.6, n = 0.43 and eV/fi = 0.002, the difference is estimated to be at most 



Figure 6 shows the eo-dependence of Gn at eV/fi = 0.002 for various values of the back- 
ground transmission T{, with {p = 0, T / fi = 0.02 and U / ^ = 0.6. At Tf, = 0.0, the conductance 
plateau is suppressed by the low bias voltage. When T5 = 0.3, the bias voltage induces a peak 
in the conductance plateau. Such a behavior is recently observed in the side-coupled QD,^ 
which is one of the systems exhibiting the Fano resonance. As the background transmission in- 
creases, the resonance peak of the differential conductance tends to split into two asymmetric 
peaks (cf. the case of T5 = 0.6). 

This could be understood as follows. Recall that the similar splitting of the resonance 
peak has been observed in the noninteracting case at higher bias voltage. In this case, a 
conductance peak is formed once the dot energy passes the Fermi level of the electrodes. At 
high bias voltage, two electrodes have very different Fermi energies and the two peaks are 
observed. While, if the bias voltage is so low that the difference between the Fermi levels 
is less than the width of the dot level, the two peaks are degenerate and one observes only 
one peak. On the other hand, as mentioned after (35), the transmission probability in the 
interacting case has the similar form as in the noninteracting case except the renormalization 
of the dot energy and the reduction of the dot level-width. Therefore, if the bias voltage 
is large enough compared with the level-width, two peaks can be observed. In the plateau 
regime, the reduction factor is small (cf. Fig. 3) and there appear two peaks even at very 
low bias voltage. In other words, because of the reduction of the width, the system becomes 
more sensitive to the bias voltage and the splitting of the resonance peak is observed at low 
bias voltage. 

The differential conductance Gn is shown as a function of the bias voltage in Fig. 7, Fig. 8 
and Fig. 9 for various occupation numbers, where the other parameters arc set again as = 0, 
T / jjL = 0.02 and U/ji = 0.6. When the background conduction is absent T), = 0.0 (i.e., in case 




1%. 
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of the transport through a single QD), the conductance has a sharp peak at zero bias voltage. 
These peaks correspond to those experimentally observed in refs.8, 9. With decreasing the 
occupation number (this corresponds to the increase of eo from the center of the conductance 
plateau of Fig.l), the zero-bias maximum is reduced from the unitary limit value 2e^//i and 
the width of the peak becomes broader (cf. the curves at n = 0.48, 0.43, 0.40). This behavior 
of the zero-bias maximum agrees with the calculations reported by Craco et al.^^ and by B. 
Dong et Note that G„ at n = 0.48 has very small negative value for larger \V\, which 

seems to be an artifact due to the use of Gn- 

In contrast, when the direct tunneling is turned on, the conductance chnages its behavior. 
At Tfc = 0.6 (Fig. 8), there appears a zero bias dip and, similarly to the zero bias peak at 
Th = 0.0, the dip width becomes broader as n decreases. Moreover, at = 0.45 (Fig. 9), a 
peak gradually changes into a broad dip as the decrease of n. In case of infinite Coulomb 
repulsion U +oo, Bulka et. al}^ predicts the possibility of the zero bias dip depending 
on the magnetic flux piercing the AB ring. Here we have shown that the similar behavior is 
possible depending on the strength of the background transmission. 

6. Conclusions 

We have investigated the transports through the QD embedded in the AB ring in the linear 
response and finite bias regimes, employing the nonequilibrium finite-?/ SBMF approach. The 
mean fields are evaluated at a nonequilibrium steady state with respect to the mean field 
Hamiltonian which is constructed with the aid of the C*-algebraic approaches. The Coulomb 
interaction is found to be taken into account as the renormalization of the dot energy and the 
reduction of the dot level-width. In particular, the latter may be interpreted as a signature of 
the formation of the singlet bound state between the dot and conduction electrons. We have 
studied the Fano-Kondo resonances and AB oscillations at the linear response regime and 
obtained the results which, by using twice larger Coulomb energy, are in good agreement with 
those by the numerical renormalization group (NRG) technique. Furthermore, this method 
is applied to the case with finite bias voltage, and we found that the Fano-Kondo plateau is 
suppressed and deformed by small bias voltage and that, in addition to the zero-bias peak, 
the zero-bias dip can be observed in the differential conductance depending on the strength 
of the background transmission. 

It is interesting to give a rough estimate of the Kondo temperature. If the SBMF param- 
eters were weakly dependent on temperature, the linear-response conductance at = and 
eo — X = fj, reads as 

for T/{z^r) ^ 1. By comparing this with the low temperature expansion of the empirical 



14/18 



J. Phys. Soc. Jpn. 



Full Paper 



formula^^ 



( 



) 



s 




T2 + (21A - i)r2 



where Tk stands for the Kondo temperature, one obtains s = 5/37 and ~ 2.64 x z'^T. In 
case of r//i = 0.02 and C///i = 0.6, ~ 0.066 and Tk ~ 375mK for /x = 9.29 meV. The value 
is consistent with Haldane's formula:^^ 



which leads to Tk — 443mK at n = 0.48 as well as the experimental value of Tk = 450 ~ 
460mK.'^ It should be noted that, without the reduction of the dot level-width, one would have 
a ten times larger value for the Kondo temperature. Moreover, this analysis shows that the 
reduction factor z'^ of the dot level-width is essentially the Kondo temperature and, thus, the 
eo-dependence of z'^ (cf. Fig. 3) can be reinterpreted as eo-dependence of the Kondo temper- 
ature, which is parabolic near the center of the plateau and is consistent with experiments.^ 
These observations seem to support our interpretation that the reduction of the dot level- 
width implies the formation of the Kondo singlet. However, to confirm this view, the finite 
temperature case should be investigated and will be discussed elsewhere. 
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Appendix: Construction of 'in' fields 

Here wc derive the 'in' fields and show (25) and (26). Thanks to the bilinearity of Hp, the 
'in' fields fikc^r can be written as 

f^kar = akar + 9^^\k)fa + X] / ^^'di^i^, k')ak'ar' ■ (A-1) 

r' 

Substituting this into [HF,(3kar] = —^krl^kar^ we obtain e.g., 

and the boundary condition: a^crri't) ex.p{iujkrt/^) har (as t — >■ — oo) gives 

gfUk.k') = (zJl\k) + We-^^ I dk"gf}^{kX)u,„^ . 

i^kL — <^k'L — \ J J 

In the same way, a set of linear equations among six functions {k), g^l,{k, k') {r, r' = L, R) 
is derived and its solution leads to 

_ , Xa(^fcr) r 

Pkcr — (^kar "r a / \ ^kr^aja 
^^aKI^kr) 

, fjll( ^I'^krUk'r ^a(,l^kr) , zlukrUk'f l^c{i^kr) \ . » o\ 

+ I dk \ 7:7—7 x«*;Vr H TiTl T^fc'af , (A-2j 

J \l^kr — '^k'r — 1^ ^(T\^kr) ^kr — '^k'f — 1^ J^cry^kr) J 

where the notations are the same as those of (25) and (26). On the other hand, because of the 
completeness of which is guaranteed by the absence of the bound states for Hp^ one has 

r 

which gives (25). Eq.(26) can be derived similarly. 
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Fig. 1. The conductance in the linear- response regime as a function of the original dot level eo/M) 
where (^^ = 0, r//i = 0.02 and the Coulomb interaction is U / = 0.6 {U /T = 30) with /i the Fermi 
energy. Solid curve is the background transmission probability T5 = 0.0, the long-dashed curve at 
Tb = 0.1, the short-dashed curve at Tf, = 0.3, the dotted curve at Ti, = 0.6 and the dash-dotted 
curve at Tb = 1.0. 




DOT LEVEL [e^ln] DOT LEVEL [e^ln] 

Fig. 2. The Coidomb interaction dependence of the conductance with (p = 0.0, Tb = 0.6. The Fano- 
Kondo plateau is formed with increasing the Coulomb interaction, the left side of the upper panel 
is for [///i = 0.02, the right side for C///X = 0.3 and the left side of the lower panel for U / ji = 0.14, 
the right side for {///x = 0.6. 
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DOT LEVEL [e^ln] DOT LEVEL [eg/n] 



Fig. 3. SBMF parameters as functions of the dot energy. The left side of the upper panel is the 
differential conductance with Tt, — 0.6, the right side is the electron occupation number and 
double occupation probability, the left side of the lower panel is the reduction factor of the dot 
level- width z"^, and the right side is the renormalized dot energy as functions of the original dot 
level. The Coulomb interaction and AB phase are set to be U/T = 0.6, cp = 0. 




Fig. 4. The conductance for various values of the AB phase: <fi = (solid curve), (fi = 7r/2 (short- 
dashed curve), ip = w (long-dashed curve). The other parameters are set to be Ti, = 0.6, U/jj. = 0.6. 
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Fig. 5. AB oscillations as a function of the AB phase for the fixed electron occupation number with 
Th = 0.3 and U/fj, = 0.6. Individual curves correspond to n = 0.2 (solid curve), n = 0.45 (short- 
dashed curve), n = 0.5 (long-dashed curve), n = 0.55 (dotted curve) and n = 0.8 (dash-dotted 
curve). 
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DOT LEVEL [Zf^n] 

Fig. 6. Differential conductance at nonequilibrium regime as a function of the original dot level for 
several background transmission probabilities Tb = 0.0, = 0.6 (left and right of upper panel) 
and U = 0.3 (left side of lower panel) with ip = 0, U/jj, = 0.6 and eV/n = 0.002. 
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Fig. 7. The differential conductance as a function of the bias voltage with Tt = 0.0, (p ~ and 
U/^ = 0.6 for several fixed occupation munbers: n = 0.48 (solid curve), n — 0.43 (long-dashed 
curve) and n = 0.40 (short-dashed curve). At n = 0.48 and eV/ ji = 0, the conductance reaches 
the unitary limit. 
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Fig. 8. The differential conductance as a function of the bias voltage with Tj, = 0.6, (fi = and 
U/jjL = 0.6 for several fixed occupation numbers: n = 0.48 (solid curve), n = 0.43 (long-dashed 
curve) and n = 0.40 (short-dashed curve). 
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Fig. 9. The differential conductance as a function of the bias voltage with Th = 0.45, if = and 
U/n = 0.6 for several fixed occupation numbers: n = 0.48 (solid curve), n = 0.43 (long-dashed 
curve) and n = 0.40 (short-dashed curve). 
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